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Abstract. We investigate applicability of the Hasselmann 
kinetic equation to the spectrum of surface gravity waves at 
different levels of nonlinearity in the system, which is mea- 
sured as average steepness. It is shown that even in the case 
of relatively high average steepness, when Phillips spectrum 
is present in the system, the spectral lines are still very nar- 
row, at least in the region of direct cascade spectrum. It al- 
lows us to state that even in the case of Phillips spectrum the 
kinetic equation can be applied to the description of the en- 
sembles of ocean waves. 



1 Introduction 

Two really seminal papers were published in the area of 
Physical Oce anography abou t half a century ag o. We mean 
the art icle by IPhillipsI (Il958h and the work of iHasselmaim 
(119621) . Both of them were concentrated on the same prob- 
lem: what is going on with ocean waves growing under in- 
fluence of wind? 

O. Phillips suggested that this growth is arrested by wave- 
breaking, in other words by formation of "white caps" (or 
"white horses"). Wave breaking is the main mechanism of 
energy dissipation. In "white caps" mechanical energy of 
waves transforms to turbulence on small scales, then to heat. 
This is a strongly nonlinear phenomenon which cannot be 
studied by perturbative methods. An analytical theory of 
wave breaking is not developed yet. The Phillips' assumption 
on predomination of wave-breaking effects made him possi- 
ble to offer, on the base of dimensional considerations, the 
universal energy spectrum of wind driven waves, so called 
Phillips' spectrum. 

K. Hasselmann developed a completely different theory. 
He noticed that a typical ensemble of ocean waves hoards 



a small parameter - average steepness (slope) of the surface 
fj,. One can define it, for instance, as = V (I Vr/p), where 7] 
is the shape of the surface. This definition has the most clear 
geometrical meaning. Another definition estimates it as av- 
erage amplitude A of the waves multiplied by wave-vector 
of the spectral peak kp. This definition is easy to use for 
experimental observations. Later we shall consider another 
definition, which involves spectral distribution function. All 
these definitions give very close values of the average steep- 
ness. Typically /Lt~ 0.1. This fact allows to assume that the 
wind-driven sea is a community of weakly interacting waves, 
which can be described statistically by the use of expansion 
in powers of average steepness /i. This is a tedious proce- 
dure because one have to expand up to order fi^. However 
K. Hasselmann coped with this hard work and derived his 
famous Hasselmann kinetic equation for squared wave am- 
plitudes. 

Then it was found that Hasselmann equation has 
exact stationary solutions (Kolmogoro v- Zakhar ov 
or KZ spectra) (ZakharoyandFilonenkol 1967allb 
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Zakharov and ZaslavskiiL Il982l: IZakharov et all Il992h 
decaying at high frequency region slower than Phillips spec- 
trum. A question arises which theory is more correct? The 
answer is given by experiment. Numerous measurements 
made in lakes and ocean showed that the real spectrum 
is a combination of weak-turbulent KZ spectra and the 
Phillips spectrum. If ujp is spectrum peak frequency, in 
the energy containing range < w < Swp KZ spectra 
are realized, while the Phillips spectrum predominates 
in the high-frequency range. This fact has the following 
explanation. Wave turbulence in the wind driven sea is a 
mixture of weak and strong turbulence. This is a question 
of phase correlation. Weakly nonlinear interaction provides 
correlation of waves phases on large scales - hundreds 
of characteristic wavelengths. At the same time the wave 
breaking and white capping is the localized phenomenon. 
As a result looking at the wind driven sea one observes 
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formation of short living localized wave breaking events 
embedded into homogeneous weakly nonlinear background. 

Similar situation, coexistence of weak turbulence and 
localized coherent structures, is typical for wave turbu- 
lence. It takes place, fo r instance, in nonlinear op- 



tics dDvachenko et all Il992h where optical turbulence co- 



exists with self-focusing or in isotropic plasma where 
weak Langmuir turbulence c oexists with Langmuir col- 



lapses (IZakharov et al.L 119921) . Similar phenomena were 



observe d recently in numeri cal simulations for gravity 



waves (IZakharov et al.L l2007h . In all these situations the 



usual course of action is to augment corresponding ki- 
netic equations of weak turbulence by introduction of an 
additional term, describing dissipation of energy in the 
coher ent structures ( corresponding to white capping on- 
sets) (IZakharov et al.l i2009). The question is it still reason- 
able to use kinetic equation as an adequate model of waves 
interaction in the presence of coherent structures? 

In this paper we give positive answer to this question. To 
do this we perform massive numerical simulation of primor- 
dial dynamical equations describing potential flow of Euler 
equations of ideal fluid with a free surface. In our numerical 
experiments we studied Fourier spectra of space-time corre- 
lation function of normal canonical variables describing the 
surface dynamics. Our goal was to determine the "shape of 
the line", i.e. frequency spectra of certain spatial Fourier har- 
monics. According to the theory of weak turbulence these 
spectra must be narrow and concentrated near the frequency 
given by the dispersion relation for waves of small ampli- 
tude. We obtained the following important result - in a broad 
range of wave numbers the spectra remain narrow even in the 
presence of wave-breaking events, when the spatial spectrum 
obeys the Phillips law. This result opens a door to construc- 
tion of a properly justified dissipation term which can be used 
thereafter in wave prediction operational models. 



2 Basic model 

We solve numerically weakly nonlinear Euler equations for 
dynamics of incompressible deep fluid with free surface 
in the pres ence of gravity by the ps eudo-spectral code de- 



scribed in Korotkevich et al. ( 20121). The code was ve r- 



ified in our previous pap ers ("Pvache nko et al. , '2003 a'b 



2004 IZakharov et all 12005. 2007; Korotkevich et al.. 2008 : 
KoroflievichL l2008ll20 1 2h . The equations are written for sur- 
face elevation T]{x,y,t) and hydrodynamic velocity potential 
on the surface ip{x,y,t). Equations are result of weakly non- 
linear expansion o f the Hamiltonian up to the fourth order 
terms in steepness dZakharov et al. , 1992 ) 
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Here dot means time-derivative, A - Laplace operator in 
r — (a;,y)-^-plane, k = A, F^^ is an inverse Fourier 
transform , 7fc i s a d issipation rate (according to recent work 
by Dias et alj (E008:) it has to be included in both equa- 
tions), which corresponds to viscosity on small scales and, if 
needed, "artificial" damping on large scales. is the driving 
term which simulates pumping on relatively large scales (for 
example, due to wind). Here and further we use symmetric 
Fourier transform 



zvr 



Normal variables are introduced as follows 



Vk - 



2LUk 



(2) 



where ojk — y/gk, g - gravity acceleration which was taken 
g = l. The equations were solved in a 27r x 2tt periodic box. 
Number of modes was 1024 x 1024. Inside the box spectra 
were practically isotropic, thus one can put 
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We define the "current steepness" fik by relation 



K 

= 27ry k'^/^Ukdk. 



(3) 



Here we use angle-averaged and summation is replaced 
by integration over k = \k\. Average steepness can be ob- 
tained as a limit: 



lim ^fe. 

k—^+oc 



This definition gives a value of average steepness which dif- 
fers from the geometrical definition fi = VvT^^F) '^^^ 
periments at max by value close to 0.005, which is just sev- 
eral percent of the characteristic average steepness. 

For the spatial harmonic with wave vector k we defined 
the time Fourier transform as follows 



a{k,uj) = i / a{k,t)e- 



*dt. 



Here T — 27r/wfe - period of the chosen wave. Thereafter 
I{k,bj) ~ |a(fc,w)p. In equations ([T]i pumping had the form 



Pk = fke 



iRhit) 



fk 



AP (fc-fcpl)(fcp2-fc) 

— if fc < kpi OTk> kp2; 



(4) 



here kpi = 28, kp2 — 32 and Fq = 1.5 x 10^^; Rk{t) was a 
uniformly distributed random number in the interval (0,27r] 
for each k and t. The initial condition was low amplitude 
noise in all harmonics. Time step was At — 6.7 x 10^**. The 
dissipation function 7^ — ^^^^ +7fe^'- 
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Fig. 1. Spectrum with both condensate and inverse cascade sup- 
pressed. Linear scale. 



Fig. 2. Spectrum with both condensate and inverse cascade sup- 
pressed. Logarithmic scale. 



Artificial viscosity 7^^^ was the same in all experiments 



0,iik<kd- 



(5) 



where kd = 256 and 70 = 0.97 x 10^. Due to presence of 
this dissipation the most part of our phase space, actually for 
k > 256, was "passive". But this was necessary in order to 
get rid of aliasing harmonics and to simulate dissipation of 

energy. 

(2) 

7^ was dissipation concentrated in small wave numbers. 
It was zero only in one of three experiments. 



3 Description of experiments 

We performed three series of experiments, choosing different 
functions 7^ ' , describing damping in the area of small wave- 
numbers. 

3.1 Without condensate and inverse cascade 

In the first series of experiments we assumed that it was lin- 
early growing in small wave-numbers fc < 28 and equal to 



(2)^ j0.2|fc-28|,fc<28, 
~|0,iffc>28. 



In this case the inverse cascade of wave action was com- 
pletely suppressed. We observed formation of the direct 
cascade, reasonably described by the standard KZ-spectrum 
jflfep ~ /3iA;~'*, with /?! ~ 1.2 x 10~*. This spectrum was ob- 
served in the range of scales 32 < fc < 150. The spectra in the 
linear and logarithmic scales are presented in Figures [T]0 
The steepness grows in this range approximately as shown 
in Figure [3] At k^ax — 200 steepness reaches its saturation 
level Umax — 0.104. The shapes of spectral lines at fc equal 
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Fig. 3. Steepness as a function of absolute vale of wave-number 
Both condensate and inverse cascade were suppressed. 



to 50, 100, 150, 200, and 250 are represented in Figures iH 
One can see that the spectral lines are narrow. Only at 
area of large k ~ 250 we observe bounded (slave) harmonics 
which do not obey hnear dispersion relation. 
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Fig. 4. Spectral line of the harmonic fe = (0,50)^. Both condensate 
and inverse cascade were suppressed. 
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Fig. 5. Spectral hne of the harmonic fe = (0,100)^. Both conden- 
sate and inverse cascade were suppressed. 
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Fig. 7. Spectral line of the harmonic fe = (0,200)^. Both conden- 
sate and inverse cascade were suppressed. 
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Fig. 8. Spectral hne of the harmonic fe = (0,250)^. Both conden- 
sate and inverse cascade were suppressed. 



Fig. 6. Spectral line of the harmonic fc = (0,150)^. Both conden- 
sate and inverse cascade were suppressed. 
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Fig. 9. Spectrum with condensate suppressed. Linear scale. 




Fig. 11. Steepness as a function of absolute vale of wave-number 
Only condensate was suppressed. 
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Fig. 10. Spectrum with condensate suppressed. Logarithmic scale. 



3.2 Without condensate, inverse cascade is present 

In the second series of experiments the low frequency damp- 
ing was presented only in small wave numbers k < 10, where 

(2) 

the dissipation rate was constant 7^ ' — 0.05. In this ex- 
periment we observed formation of both direct and inverse 
cascades. The dynamic range for inverse cascade was too 
short for accurate evaluation of the power of the spectrum, 
the only statement which can be made is that it is relatively 
close to the predicted \ak\'^ ^ k^'^'^/^ spectrum, with power 
exponent a in A:^" ranging from 3.1 to 4.0. In the area of 
direct cascade we observed slightly steeper spectrum than 
KZ-spectrum slope. In the area 32 < A: < 150 the spectrum 
can be approximated by |afcp ^ k~'^-^. This deviation from 
pure KZ-slope /c"^^ can be explained by the influence of wave 
breaking and white capping effects (Korotkevich, 2 012) . The 
spectra are presented in Figures |9][T0] The steepness in 
this case reaches the level /i~ 0.130 (see Figure fTTll. The 
spectral lines in this series of experiments were still narrow 



Fig. 12. Spectral line of the harmonic k — (0,50)^. Only conden- 
sate was suppressed. 



for k — 50,100,150,200. In the area of significant damping 
k — 250 we observed intensive formation of slave (bond) har- 
monics (see Figures n"2lfT6l l. 
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Fig. 13. Spectral line of the harmonic fc = (0, 100)^ . Only conden- 
sate was suppressed. 
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Fig. 14. Spectral hne of the harmonic fc = (0, 150)^. Only conden- 
sate was suppressed. 



Fig. 16. Spectral line of the harmonic fc = (0,250)^. Only conden- 
sate was suppressed. 
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Fig. 15. Spectral line of the harmonic fe = (0,200)^. Only conden- 
sate was suppressed. 
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Fig. 17. Spectrum with condensate. Two grid sizes. Linear scale. 
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Fig. 18. Spectrum with condensate. Two grid sizes. Logarithmic 
scale. 



3.3 With both condensate and inverse cascade 

In the last series of experiments we completely eliminated 

(2) 

low-frequency dissipation 7j, — 0. In this case we observed 
formation of intensive inverse cascade leading to creation of 
the condensate at fc ~ 5 — 4. The spectrum in the area of in- 
verse cascade was the same as in the previous series of exper- 
iments, while in the area of direct cascade we observed the 
Phillips spectrum lofeP ~ ^7^^^ i nstead of a weak-turbulent 
KZ-spectrum ( K orotkevichl 20081) . The spectra are presented 
in Figures [T7][T8] In order to check that spectrum slope 
is not changing any more, long time calculation was per- 
formed on a smaller grid, which showed motion of the con- 
densate position by one wave-number grid step without any 
significant influence on the inverse cascade region. The av- 
erage steepness was essentially higher than in the previous 
experiments reaching its "limiting value" /i ~ 0.142 (see Fig- 
ure [T9]l. We would hke to stress that the difference with 



Fig. 19. Steepness as a function of absolute vale of wave-number. 
Condensate was present. 



the previous case, which looks quite small in absolute value 
(0.130 and 0.142), is really quite strong, because the dissi- 
pation rate depe nds strongly on the ave rage steepness in this 



range of values (IZakharov et al.L l2009h as well as probabil 



2000). 



ity of white capping grows pretty fast (IBanner et al 
It is important to notice that the steepness of the condensate 
was quite moderate /i ~ 0.06 (it is worth to note, that av- 
erage steepness getting significant contribution from small 
scales, quite far from the spectral peak, which means that of- 
ten used definition of the average steepness through the prod- 
uct of mean amplitude and wave-number of spectral peak can 
deviate significantly from the geometrical definition through 
the average slope of the surface). However modulation of the 
short waves by long waves caused intensive micro-breaking 
of waves which forms t he Phillips spectrum (the process 
described qualitatively in iKorotkevichI (2008, 2012h ). The 
most interesting and important question is about the shapes 
of spectral lines in the area of Phillips spectrum. They are 
presented in Figures l20ll24l For this experiment we analysed 
the longest time series which resulted in higher frequency 
resolution. After obtaining the frequency spectrum we used 
moving averaging in order to get rid of noise. One can 
see that in the most interesting area 32 < A: < 150 the spec- 
tral lines are still narrow while essentially "contaminated" 
by the slave harmonics. The maxima of the spectral peaks 
are shifted to the high-frequency area according to the theo- 
retically predicted nonlinear frequency shift. In the area of 
shorter waves k > 200 the spectrum is a chaotic mixture of 
leading and slave harmonics. 



4 Conclusions 

We believe that our experiments make possible a "marriage 
of Heaven and Hell" in spirit of William Blake. Both out- 
standing scientists - O. Phillips and K. Hasselmann are right. 
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Fig. 20. Spectral line of the hannonic k = (0,50)^. Condensate 
was present. 



Fig. 22. Spectral line of the harmonic fc = (0,150)^. Condensate 
was present. 
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Fig. 21. Spectral line of the harmonic k = (0,100) . Condensate 
was present. 



Fig. 23. Spectral line of the harmonic fe = (0,200)^. Condensate 
was present. 



If the local steepness /U is small (/z < 0.1) the Hasselmann 
kinetic equation is valid without any augmentation by any 
additional dissipation terms. For /z significantly higher than 
0.1, in the area of very short waves the kinetic equation is 
not appUcable in its "pure" form. It is impossible to separate 
"leading" and "slave" harmonics in this area. This part of the 
ocean spectrum cannot be described analytically by the use 
of perturbative methods. Nevertheless, the Phillips spectrum 
in this area is still applicable. The only theoretical reason 
for this statement is dimensional considerations supported by 
experimental data. At the same time we claim that there is 
a "grey area", which is especially interesting, because it is 
containing most of energy and generate most of steepness, 
where micro-breaking have equal foot with weakly nonUn- 
ear resonant waves interaction. Spectra in this area could be 
described by an "augmented" Hasselmann equation, includ- 
ing an additional term describing dissipation energy due to 



wave-breaking. Similar additional "dissipation terms" S^iss 
are widely used in well developed operational wave forecast- 
ing models. But they are introduced "out of the blue" and not 
supported neither by theoretical consideration nor by exper- 
imental observations. Moreover, it is not clear a priori that 
one can use the Hasselmann kinetic equation in a situation 
where wave-breaking events are frequent enough. We hope 
that our experiments showed that this is possible. A nec- 
essary and sufficient condition of applicability of the Has- 
selmann equation is narrowness of the spectral line. In the 
present article we assert that in the "grey area" with a fre- 
quency range in one half of the decade the frequency spectra 
of harmonics are still narrow lines. It means that the Has- 
selmann equation is applicable there. Of course it must be 
augmented by a proper dissipative term. The existing and 
widely used dissipative terms hardly can be correct. They 
do not satisfy to the minimal requirement - the Hasselmann 
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Fig. 24. Spectral line of the harmonic k — (0,250)"^. Condensate 
was present. 



equation in the presence of this term has to have the Phillips 
spectrum as an exact solution. Hence, the urgent problem 
is a construction of a "really good" dissipative term. One 
po ssible and hopefully pla usible variant was recently offered 
by IZakharov et al.l (120121) . The dependence of the dissipa- 
tion term on the average steepness was recently measured di- 
rectly from the n umerical experiment, th e preliminary results 



can be found in IZakharov et al. I (|2009'). In our next paper 



we shall analyse what information about possible shape of 
the dissipative term can be extracted from massive numerical 
experiments. 
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